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We investigate the electronic structure of the frustrated magnet Na^^Os using density functional 
theory. Due to strong spin-orbit coupling, the hyperkagome lattice is characterized by a half-filled 
complex of d states, making it a cubic iridium analogue of the high temperature superconducting 
cuprates. The implications of our results for this unique material are discussed. 

PACS numbers: 71.70.Ej, 75.10.Lp, 75.30.Et, 75.50.Ee 



I. INTRODUCTION 

Recently, several materials have been identified^ as 
possible candidates for a long sought state of matter - the 
quantum spin liquid. - One of these candidates, Na^raOs, 
is characterized by a three dimensional cubic space group, 
where the iridium sites form a hyperkagome lattice of 
corner sharing triangles. This remarkable material is an 
insulator with a large Curie- Weiss temperature (650 K) 
and a large effective moment (1.96 /^b), yet does not 
exhibit any sign of magnetic order down to the lowest 
measured temperature^ 

There have been a number of theoretical studies of 
this iridate^ which attempt to address the nature of the 
spin liquid ground state. These studies invariably assume 
that the material is described by an isotropic Heisenberg 
model, since anisotropy acts to stabilize magnetic order. 
Isotropy would be a surprise given the strong spin-orbit 
coupling for the iridium ions. But in an exhaustive study, 
Chen and Balents^ formulated under what conditions 
such an effective model might be realized. To address 
various issues raised in that paper, we calculated the elec- 
tronic structure of this material within the local density 
approximation, then performed a tight binding fit in or- 
der to estimate the various exchange integrals which enter 
the spin Hamiltonian for this frustrated magnet. We find 
that the spin Hamiltonian should be strongly anisotropic, 
but offer two possible scenarios where isotropy might be 
restored. 



II. CALCULATIONAL DETAILS 

Although Na 4 Ir 3 08 has a cubic space group, the 
unit cell is quite complicated, comprising four formula 
units (60 atoms).— Moreover, the lattice breaks inver- 
sion symmetry. The space group, P4!32, is also found in 
the non-centrosymmetric superconductors L^PdsB and 
Li2Pt3B£ The material is formed from distorted IrOg 
octahedra, each comprised of two different oxygen sites - 
four of type 2 (02), two of type 1 (01). The latter should 
not be thought of as 'apical' oxygens in the cuprate sense, 
since they are not related by inversion through the irid- 
ium site. Rather, each octahedron is described by a C2 
axis along a (110) direction^ There are twelve of these 



TABLE I: Coordination shells around an iridium ion in the 
Na4lr308 structure out to a distance of 0.49a. Listed is the 
type of ion, their number, (approximate) direction, and dis- 
tance (in units of the lattice constant, a). 



type 


number 


direction 


distance 


02 


2 


100 


0.2273 


02 


2 


100 


0.2280 


01 


2 


100 


0.2333 


Ir 


4 


110 


0.3464 


Na3 


2 


110 


0.3501 


Na3 


2 


110 


0.3587 


Nal 


2 


110 


0.3612 


Na3 


1 


110 


0.3627 


Na2 


1 


110 


0.3683 


02 


2 


111 


0.4065 


02 


2 


111 


0.4225 


01 


2 


111 


0.4375 


02 


2 


111 


0.4838 



axes, one associated with each of the twelve iridium ions 
in the unit cell. Each oxygen ion is also at the center of 
a distorted octahedron, with 01 octahedra of the form 
OIr3Na3, and 02 octahedra of the form OIr2Na4. Each 
oxygen ion is also surrounded by a very distorted ar- 
ray of twelve other oxygen ions along approximate (110) 
directions. Finally, each iridium ion is surrounded by 
four other iridium ions, also along approximate (110) di- 
rections. These form the corner sharing triangles that 
characterize the hyperkagome lattice. A list of the first 
few coordination shells around an iridium ion is given in 
Table I. 

When considering an electronic structure calculation 
for this material, the first issue is that there are three 
different types of sodium ions, two of which are partially 
occupied with a stochiometry of 75%£ Therefore, to form 
a simple unit cell, some approximation has to be made for 
these other two sodium sites. We choose to fully occupy 
the 12d sites (sodium type 3) and replace the 4a sites 
(sodium type 2) by empty spheres. The resulting lattice 
has 64 sites in the unit cell - 12 Ir, 8 01, 24 02, 4 Nal, 
4 empty (i.e., Na2), and 12 Na3. 

The calculations in this paper were performed using 
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FIG. 1: (Color online) Density of states (DOS) in states/Ryd 
for the scalar relativistic calculation of NaziIrsOs. The Fermi 
energy, Ef, is at zero. 



the local density approximation (LDA) within a linear 
muffin tin orbital schemed Overlap (combined) correc- 
tions were incorporated, but the Coulomb potential at 
each site was treated in a spherical approximation. The 
exchange-correlation potential used was that of Hedin 
and Lundqvist. 5 Noting that the cubic lattice constant 
is 8.985A, a muffin tin radius of 1.238A was used on the 
iridium sites, and 1.425A was used on all other sites. 
These radii were chosen based on the various near neigh- 
bor separations. s,p basis functions were used on all sites 
(including the empty ones) except for the iridium sites, 
where d basis functions were also included. For the scalar 
relativistic calculation (all relativistic effects but spin- 
orbit), this resulted in a secular matrix of dimension 316, 
which was double this (632) for the full relativistic (Dirac 
equation) calculation. The calculations were converged 
using 56 k points in the irreducible wedge of the Brillouin 
zone. We should note that the Brillouin zone for this ma- 
terial is simple cubic. Though inversion symmetry is bro- 
ken, one can restrict the irreducible wedge to l/48th of 
the zone. The eigenvalues at a general k point are doubly 
(Kramers) degenerate unless spin-orbit is incorporated, 
where this degeneracy is broken. After convergence, 165 
k points were generated in the irreducible wedge. The 
eigenvalues were then fit using a Fourier spline series^ 
with 560 lattice periodic functions. A linear tetrahedron 
methodic was used to evaluate the density of states, with 
the zone broken down into 48 ■ 8 5 tetrahedra. 



III. RESULTS 

We first present results from the scalar relativistic cal- 
culation, with the density of states shown in Fig. 1. The 
01 states form a narrow complex of bands ranging from 
-8.2 to -6.7 eV, followed by the wider 02 complex ranging 
from -6.7 to -3 eV (there is only a narrow gap between 
these two complexes). There is a gap of 1.9 eV to the t2 ff 
complex of iridium d bands ranging from -1.1 eV to 0.4 
eV. This is followed by a crystal field gap of 1.9 eV to an 
e g complex ranging from 2.3 eV to 3.6 eV. For energies 
above this, the states are primarily from the sodium sites. 
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FIG. 2: Energy bands of the t2 g complex for the scalar rela- 
tivistic calculation of Na 4 Ir 3 8 . T = (0,0,0), X = (1,0,0), 
M — (1, 1, 0) and R — (1, 1, 1) in units of ir/a. The horizontal 
line marks Ef- 
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FIG. 3: (Color online) Decomposition of the (a) t2 3 states 
and (b) e 9 states, averaged over each band, into their individ- 
ual components relative to the C2 axis of the distorted IrOe 
octahedra (scalar relativistic calculation). 

The main complex of interest is the t2 9 one (Fig. 2). It is 
split into an array of 28 bands (each Kramers degenerate) 
separated by a gap of 0.2 eV to an upper array of 8 bands. 
This basic splitting of 28 and 8 arises from the t p d hop- 
ping integral, as identified from tight binding fits, and is 
related to the four-fold degeneracy of the eigenvalues at 
the R point of the zone. There is a crystal field splitting 
of the t2 g states which was identified by projecting onto 
the three crystal field components due to the C2 axis, 
with the xy states (for a 110 axis) lying lowest in energy, 
but it is rather ill defined (Fig. 3) and is small relative to 
the overall t2 g bandwidth. The t2 S complex is composed 
of 70% Ir 5d character, 20% 02 2p character, and 4% 01 
2p character. The ccntroid of the e g complex relative to 
the t2 g one is 3.5 eV, with the t2 g -e g admixture (due to 
the low iridium site symmetry) quite small. The density 
of states at the Fermi energy (per Ir per spin) is 3.35 
states/eV. The Fermi surface (not shown) is composed 
of small pockets, two about the R point (electron- like), 
two about the M point (hole-like), and one around the 
r point (hole- like). 

Spin-orbit coupling has a profound impact on these re- 
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FIG. 4: (Color online) Density of states (DOS) in states/Ryd 
for the full relativistic calculation of Na4lra08. 
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FIG. 5: Energy bands of the t2 9 complex for the full rela- 
tivistic calculation of Na^raOg. The horizontal line marks 
Ef ■ 

suits, as can be seen in the remaining figures. The t2 g 
states now split (Figs. 4 and 5) into a lower complex of 
48 bands and an upper complex of 24 bands (the num- 
ber of bands are now doubled because of the coupling of 
inversion breaking with spin-orbit), with the e g complex 
composed of another 48 bands. The upper ta g complex is 
half- filled, with the Ir d states mostly of j — 5/2 charac- 
ter, with weak admixture of j = 3/2 character (Fig. 4). 
As explained by Chen and Balents^ and found as well 
for the electronic structure of Sr2lr04fii this behavior is 
due to crystal field admixture of the j = 5/2 and j = 3/2 
states, which acts to form a mixed lower (t2 9 ) quartet 
of states (i.e., four per iridium site), a mixed (e g ) upper 
quartet, with an unmixed (t2 S , j — 5/2) doublet in be- 
tween that is half-filled. This doublet is predicted to have 
83% j = 5/2, raj = ±3/2 and 17% j = 5/2, rrij = ±5/2 
character, which is in good agreement with the results 
we find for the lowest 8 bands of this complex (Fig. 6). 
The upper 16 bands have a slightly reduced rrij = ±3/2 
weight due to weak admixture with rrij = ±1/2 states as 
well as j = 3/2 states, which is a consequence of the low 
iridium site symmetry. Again, the density of states at 
the Fermi energy is not large - per Ir per spin it is 1.51 
states/eV. The Fermi surface is composed of three elec- 
tron like sheets around T, and two hole-like sheets around 
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FIG. 6: (Color online) Decomposition of the states of the up- 
per %2 g 'doublet' complex, averaged over each band, into their 
atomic spin-orbit components, j, ±rn.j (full relativistic calcu- 
lation). The weights are normalized to the total d weight. 
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FIG. 7: Fermi surface for the full relativistic calculation of 
Na4lr308 in the symmetry planes of the zone. 



the R point (Fig. 7). The bands are doubly degenerate 
on the zone face due to the non-symmorphic nature of 
the space group. 

We now wish to connect these results to the insulating 
behavior of the actual material. In the case of Sr2lr04, 
inclusion of the Coulomb repulsion is not sufficient to 
open an energy gap unless spin-orbit is included. In 
essence, spin-orbit produces a half-filled doublet which 
can then 'Mott-Hubbardize' and form an insulator, sim- 
ilar to what occurs in the cupratesJ^ Only a small U is 
needed to cause thisJaL Similar considerations should ap- 
ply to Na4lr308- From the centroids of the two quartets 
and the doublet, we have inferred the parameters Acf 
(splitting of t2g and e 9 ) and A (coefficient of I ■ s for the 
spin-orbit coupling) of 3.43 eV and 0.58 eV, respectively. 
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IV. TIGHT BINDING ANALYSIS 

To investigate these issues further, we performed tight 
binding fits to the electronic structure, restricting to Ir 
5d and O 2p orbitals. This involved generating Slater- 
Koster parameters,— and then performing a least squares 
fit to the band eigenvalues using Powell's method^ Be- 
cause of the distorted nature of the lattice, there are a 
large number of tight binding parameters, even if one re- 
stricts to near neighbors. This is because there is not 
a unique distance for several of the atom combinations 
(Ir-02, 01-01, 01-02, and 02-02). To compensate for 
this, we assumed that the tight binding parameters for 
each atom combination vary as the fourth power of the 
distance. This assumption works well for the 0-0 hop- 
pings in IrOsr^ (note from Table I that the two Ir-02 
distances are almost identical). This reduced us to four 
t p d parameters (a and tt for Ir-01 and Ir-02), three tdd 
parameters (a, tt, and S for Ir-Ir), six t pp parameters (<r 
and tt for 01-01, 01-02, and 02-02), four on-site ener- 
gies (2p energies for 01 and 02, and t 2g and e g energies), 
and the spin-orbit coupling of the Ir d (A), for a total of 
18 parameters. We have explored putting in the residual 
crystal field splittings due to the low Ir site symmetry 
(Fig. 3), but this did not lead to improved fits. 

We first converged the fit for the scalar relativistic cal- 
culation, using 42 eigenvalues (top and bottom of the 01 
complex, 02 complex, e g complex, and all 36 t 2g bands) 
at the four symmetry points (r, X, M, and R). As input, 
we used Ir02 hopping parameters^ scaled to the sodium 
iridate distances, as well as on-site energies estimated 
from averages extracted from the I decomposed density 
of states. The fit (not shown) yielded a band structure 
in good agreement with Fig. 2. We then used these fit 
parameters as input to fitting the full relativistic calcu- 
lation, with the same set of eigenvalues (now numbering 
78) at each of the 4 k points. The fit gives a good descrip- 
tion of the bands and Fermi surface, as can be seen in 
Figs. 8 and 9. These fit values are listed in Table II. One 
can see that the t p d and tdd parameters have comparable 
magnitude, and that the 01 and 02 integrals are quite 
different. Note that the difference of the t 2g -e g energies 
is only 2.60 eV. The difference from the centroid splitting 
of 3.43 eV quoted earlier is due to the contribution of the 
hoppings to the splitting. 



V. EXCHANGE INTEGRALS 

As we are in the strong spin-orbit coupling limit in the 
sense of Chen and Balents^ we can use the expressions 
in their paper to estimate the various exchange integrals 
(Table III). We caution that these values are very sensi- 
tive to the tight binding parameters, as well as the values 
of the Coulomb U. The exchange Hamiltonian is of the 
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FIG. 8: Energy bands for the upper t2 S 'doublet' complex 
from the tight binding fit to the full relativistic calculation of 
Na^^Os. The horizontal line marks Ef- 
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FIG. 9: Fermi surface from the tight binding fit to the full 
relativistic calculation of Na4lr308- 

form 5 

Hex = \Jd ~\~ Js)S n ■ S m + D • (S n x S m ) + S n • r • S m 

(1) 

where Jd is the direct exchange between iridium ions, J s 
the isotropic component of the superexchange mediated 
by the oxygen ions, D the Dzyaloshinski-Moriya interac- 



TABLE II: Tight binding hopping parameters in eV. The on- 
site energies are eoi = -5.1913, eo2 = -4.0954, e t2g = -1.7983, 
e e = 0.7987, with the spin-orbit coupling A = 0.6386. 





a 


7T 


5 


Ir-Ol 


-2.5579 


0.3186 




Ir-02 


-2.1070 


1.1817 




Ir-Ir 


-0.6372 


0.0719 


0.1545 


Ol-Ol 


0.6156 


0.0354 




02-02 


0.5463 


-0.2121 




01-02 


-1.1327 


0.0663 
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TABLE III: Exchange integrals in meV. The spin Hamiltonian 
is given in Eq. 1. Quoted are values where Ir site m is along an 
(0,1,-1) direction relative to Ir site n. D are the Dzyaloshinski- 
Moriya, and Y the anisotropic superexchange terms (in the 
second row for F, the jk refer to the parenthesis, xy, etc.). In 
addition, the direct exchange is Jd — 20.1, and the isotropic 
superexchange term is J a = 12.9. The last row, Ju, is the 
total exchange for the diagonal components (Jd + J s -+- r«). 
The assumed value of U is 0.5 eV— (with U p set to zero). 



i (jk) 


x (xy) 


V (xz) 


z (yz) 


D t 


47.3 


1.3 


-4.7 


r H 


36.1 


-36.8 


-36.2 


Ij* 


2.1 


-8.3 


-0.2 


J% i 


69.2 


-3.7 


-3.1 



tion, and V the anisotropic part of the superexchange. 

Remarkably, for the direct Ir-Ir exchange, we find that 
even for the distorted lattice, the interaction is isotropic, 
and is given by the simple formula Jd — 2t 2 /U, where 
* = \^dd + \^kd + ~n^dd- The value we find is sizable, 
about 20 meV, which is not so different from the exper- 
imental estimate of 28 meV4 It would be larger, except 
for the anomalously large value of t s dd . In that context, 
we remark that these hoppings are effective values, and 
compensate for the fact that we only include 18 tight 
binding parameters in our fits. For instance, in the more 
complete analysis Mattheiss did for the simpler L-O2, t dd 
is only about 6% of t dd ~ We have done a tight binding 
fit including the residual crystal field splittings (Fig. 3) 
which did have a much reduced value for t dd (with a 
larger value for Jd of 37 meV), but the Fermi surface was 
degraded in this fit. 

We now turn to the superexchange, which is respon- 
sible for all values of the exchange Hamiltonian except 
for Jd- The values (Table III) are large and highly 
anisotropic (we should caution, though, that these values 
will be further affected by the t pp hoppings 16 ). This is 
unlike what would occur if the local site symmetry of the 
Ir ions was purely cubic <H The large values are not only 
due to the distortion of the octahedra, but also because 
the superexchange between two given Ir ions is mediated 
by an 01 and an 02 ion, each of which is in a different 
environment (Fig. 10). For instance, the non-zero values 
of J,, D, and the off-diagonal components of T are due 
to the distortion, with the magnitude strongly influenced 
by the difference between the 01 and 02 ions. On the 
other hand, the reversed sign of the diagonal components 
of T relative to the cubic model of Chen and Balents^ is 
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due to the the difference between 01 and 02, with the 
magnitude little affected by the distortion. 

There are only two ways we know to reduce these large 
values of the superexchange. First, in the limit of large 
Up (Coulomb repulsion on the oxygen sites), they would 
vanish. Second, because the sodium sites about the oxy- 
gens are only partially occupied in the real material, it 
is possible that the resulting randomness would dephase 




FIG. 10: (Color online) Superexchange pathway between two 
Ir ions as marked by the arrows. The Ir ions are the large 
(red) spheres and the oxygen ions the small ones. The lower 
left (green) pathway is via an 02 (blue) ion, the upper right 
(purple) pathway via an 01 (brown) ion. 



the various superexchange pathways. Unless one of these 
mechanisms is effective, the total exchange is predicted 
to be anisotropic. This would cast doubts on theories 
for Na^^Os based on an isotropic Heisenberg model. 
Regardless, we believe that the tight binding parameters 
we provide here will be useful for future modeling of this 
fascinating material. In particular, given the resemblance 
of its electronic structure to that of cuprates, one might 
wonder whether a doped version of this material would 
be a high temperature superconductor. 
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